A general variational principle for spherically symmetric perturbations in 

diffeomorphism covariant theories 



O 

o 



Michael D. SeifertQ and Robert M. Walcfl 
Enrico Fermi Institute and Department of Physics, 
University of Chicago, 5640 S. Ellis Ave., Chicago, IL, 60637, 



USA 



We present a general method for the analysis of the stability of static, spherically symmetric so- 
lutions to spherically symmetric perturbations in an arbitrary diffeomorphism covariant Lagrangian 
field theory. Our method involves fixing the gauge and solving the linearized gravitational field 
equations to eliminate the metric perturbation variable in terms of the matter variables. In a wide 
class of cases — which include f(R) gravity, the Einstein-aether theory of Jacobson and Mattingly, 
and Bekenstein's TeVeS theory — the remaining perturbation equations for the matter fields are sec- 
ond order in time. We show how the symplectic current arising from the original Lagrangian gives 
rise to a symmetric bilinear form on the variables of the reduced theory. If this bilinear form is 
positive definite, it provides an inner product that puts the equations of motion of the reduced 
theory into a self-adjoint form. A variational principle can then be written down immediately, from 
which stability can be tested readily. We illustrate our method in the case of Einstein's equation 
with perfect fluid matter, thereby re-deriving, in a systematic manner, Chandrasekhar's variational 
principle for radial oscillations of spherically symmetric stars. In a subsequent paper, we will apply 
our analysis to f(R) gravity, the Einstein-aether theory, and Bekenstein's TeVeS theory. 
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I. INTRODUCTION 

In recent years, there have been a number of propos- 
als for new theories which modify general relativity in an 
attempt to either explain outstanding experimental prob- 
lems or provide a theoretical framework for the analysis 
of new phenomena. In the first category, we have theories 
such as Carroll et al.'s f(R) gravity [lj, which attempts 
to modify general relativity to explain the cosmic accel- 
eration, and Bekenstein's TeVeS theory Q, an attempt 
to formulate a covariant version of Milgrom's MOND. In 
the second category, we have theories such as Jacobson 
and Mattingly 's "Einstein-aether theory" 0, a toy theory 
in which Lorentz symmetry is dynamically broken by a 
vector field which is constrained to be unit and timelike. 

If such theories are to be phenomenologically viable, 
they must satisfy two criteria, among others. First, they 
must possess solutions which are quasi-Newtonian, i.e., 
they must reproduce the dynamics of Newtonian gravity 
up to some small relativistic corrections. Second, these 
quasi-Newtonian solutions must be stable, or at least 
must not be unstable on time scales sufficiently short to 
interfere with known physical phenomena in Newtonian 
gravity. In particular, there must exist static, spherically 
symmetric solutions of the theory, corresponding to the 
interior and exterior of a star, and moreover these solu- 
tions must not be unstable on a time scale significantly 
shorter than cosmological time scales. 

It would therefore be useful to obtain a general method 
by which the stability of an arbitrary Lagrangian field 
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theory might be analyzed. The perturbational equations 
of motion off of some background will be a set of linear 
partial differential equations in terms of some perturba- 
tional fields. We shall denote these fields tp a , where a 
is an index running over the collection of perturbational 
fields. Suppose that it were the case that these equations 
take the form 



dt 2 







(1) 



where T, the time-evolution operator, is a linear opera- 
tor containing spatial derivatives only. Suppose, further, 
that we can find an inner product (•, •) on the vector space 
of fields tp a in which the operator T is self-adjoint, i.e., 
the domain of T coincides with the domain of and 



(<P, ?X) = (Tip, X) 



(2) 



Let ujq denote the 



for all ip and x m the domain of T 
greatest lower bound of the spectrum of T. Then it can 
be shown [J] that if ujq > then the background solution 
is stable, whereas if Uq < then perturbations exist that 
grow exponentially on a timescale r = l/|wo|. Further- 
more, for all if) in the space of perturbational fields, we 
have (by the Rayleigh-Ritz principle) 



< (T) % 



(jW) 



(3) 
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Thus, in practice, then, we can obtain a bound on the 
spectrum of T by plugging various "trial functions" , ip a , 
into ([3]) . If we find a trial function which yields a negative 
result, we know that there exist solutions of ([TJ) that grow 
exponentially in time, and we also obtain a upper bound 
on the the timescale on which the instability of fastest 
growth occurs. 
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Unfortunately, for a generally covariant theory, the 
perturbation equations can never directly arise in the 
simple form of ([I]). There will be gauge freedom, so the 
equations will not even be deterministic. There will be 
constraint equations, which typically will be of lower or- 
der in time derivatives than the other equations. Further- 
more, even if one succeeds in suitably fixing a gauge and 
solving some subset of the equations so that the remain- 
ing equations take the form of (pj, there may not exist 
an inner product that makes T sclf-adjoint, and even if 
such an inner product exists, it may be highly nontrivial 
to find it. 

In this paper, we will show that all of the above prob- 
lems can be solved in a wide class of diffeomorphism co- 
variant theories in the context of spherically symmetric 
perturbations of static, spherically symmetric solutions. 
The easiest task is to fix the gauge for such perturbations, 
and we make a choice of gauge in Section Mil In Section 
IIV1 we analyze the linearized constraint equations. For 
spherically symmetric perturbations, there are two inde- 
pendent constraint equations, namely a "time-time" and 
a "time-radial" equation. One of the main results of this 
paper is our proof that these two equations can always 
be reduced to a single equation of lower differential or- 
der. Furthermore, if the original constraint equations are 
partial differential equations of second (or lower) differ- 
ential order, we show that in a wide class of theories this 
resulting single equation can be solved algebraically for 
one of the metric perturbation variables, thus completely 
eliminating the constraints. In addition, we show that in 
a wide class of theories, the remaining independent com- 
ponent of the gravitational field equations can be used 
to eliminate the remaining metric perturbation variable. 
We thereby reduce the theory to solving a system of par- 
tial differential equations for just the matter perturbation 
variables. 

In the theories that we shall consider here and in a 
subsequent paper, the equations of motion of this re- 
duced theory take the form of (Q]). However, we still 
have the problem of determining whether there exists 
an inner product that makes T self-adjoint, and find- 
ing this inner product if it does exist. In this regard, 
it might be expected that a Hamiltonian formulation of 
the theory would be useful, since for a Hamiltonian of a 
simple "kinetic plus potential" form, the kinetic energy 
expression will provide such an inner product. Now, by 
assumption, the original theory was obtained from a La- 
grangian, and, from this Lagrangian, one can obtain a 
Hamiltonian [5(. Prior to gauge fixing and solving the 
constraints, a Lagrangian and Hamiltonian for pertur- 
bations can be obtained by expanding about the back- 
ground solution and keeping the terms quadratic in the 
perturbation variables. However, it is well known that 
if one substitutes a choice of gauge into a Lagrangian or 
Hamiltonian, it no longer functions as a Lagrangian or 
Hamiltonian. Similarly, if one substitutes a solution to 
some of the field equations (such as the constraints) into 
a Lagrangian or Hamiltonian, it no longer functions as a 



Lagrangian or Hamiltonian. Thus, it might appear that 
the existence of a Lagrangian or Hamiltonian formula- 
tion of the original theory will be useless for obtaining 
a Lagrangian or Hamiltonian formulation of the reduced 
theory, and also useless for finding an inner product in 
the reduced theory that makes T self-adjoint. 

In Section IVI1 we make use of the fact that — although 
substitution of a gauge choice and/or a solution to 
some of the field equations into a Lagrangian or Hamil- 
tonian does not yield a Lagrangian or Hamiltonian — 
substitution of a gauge choice and/or a solution to some 
of the field equations into the symplectic current still 
yields a symplectic current that is conserved when the 
remaining linearized equations of motion are satisfied. 
We show that — in a wide class of cases — the resulting 
symplectic current of the reduced theory gives rise to 
a bilinear form that makes T symmetric. If this bilin- 
ear form is positive definite, it provides the desired inner 
product, thus enabling us to obtain a variational princi- 
ple for determining stability. 

In Section IVIII we shall illustrate our method by ap- 
plying it to the case of the Einstein-fluid system. Wc 
thereby will re-derive in a systematic fashion the varia- 
tional principle for analyzing radial oscillations of spher- 
ically symmetric stars in general relativity that was first 
obtained by Chandrasekhar 6,7]. In a subsequent paper, 
the method will be applied to analyze f(R) gravity, the 
Einstein-aether theory of Jacobson and Mattingly, and to 
Bekenstein's TeVeS theory. 

Our notation and conventions will generally follow [§[ . 
Where relevant, we will use units in which c = G — 
1. Boldface symbols will be used to denote differen- 
tial forms. For definiteness, we will restrict considera- 
tion here to four-dimensional spacetimes that are spheri- 
cally symmetric in the sense of having an SO (3) isometry 
subgroup whose orbits are two-spheres. However, all of 
our analysis and results generalize straightforwardly to 
n-dimensional spacetimes that are spherically symmetric 
in the sense of having an SO(n — 1) isometry subgroup 
whose orbits are (n — 2)-spheres. 

II. LAGRANGIAN FORMALISM 

In this section, we shall briefly review some basic defi- 
nitions and constructions in theories derived from a dif- 
feomorphism covariant Lagrangian. We refer the reader 
to [9( and [5| for further details and discussion. 

Consider a diffeomorphism covariant Lagrangian four- 
form C, constructed from dynamical fields W that con- 
sist of the spacetime metric and perhaps additional mat- 
ter fields, assumed to be described by tensor fields. For 
computational convenience, in this paper we will use the 
inverse metric g ab rather than the metric g a b as the inde- 
pendent dynamical variable describing the gravitational 
degrees of freedom. The first variation of C can be writ- 
ten in the form 

SC = + d9 , (4) 
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which defines not only the Euler-Lagrange equations of 
motion, £ = 0, but also the symplectic potential three- 
form, 0(^f,S^). In |4]), any tensor indices of 6*f? are un- 
derstood to be contracted with the corresponding dual 
tensor indices of £, and, of course, the sum over all fields 
comprising ^ is understood. The antisymmctrized vari- 
ation of the symplectic potential yields the symplectic 
current three-form (jj(^f; Si 1 ^, <5i^), defined by 

u) = Sid(^,s 2 ^)-s 2 e(^,Si^). (5) 

(Here, 6± and S 2 are taken to be variations along a two- 
parameter family, and hence they commute, SiS 2 — 6 2 Si.) 
Applying d to this equation and using ((3]), we obtain 

do; = -(Sx£)5 2 ^+ {5 2 £)5 1 ^ . (6) 

Thus, the symplectic current is conserved, da; = 0, when- 
ever and S 2 ^/ satisfy the linearized equations of mo- 
tion. 

The Noether current three-form associated with an 
arbitrary vector field £° is defined by 

J e = 0(9, (7) 

where the "•" denotes the contraction of £ a into the first 
index of C Applying d to this equation, we obtain 

dJ 6 = -££^. (8) 

Thus, J is conserved, dJ = 0, whenever \& satisfies the 
equations of motion (i.e., £ = 0) or whenever £ a is a 
symmetry of \t (i.e., = 0). Finally, we note that a 
simple calculation Q shows that 

S3 6 = ■£ + u('$: 5-$, +d{£-0). (9) 



For the static background solution, $ and A can be 
chosen to be independent of t, in which case the only 
remaining gauge freedom is t — ► ct for some constant c. 
Thus, for the background solution, $ is unique up to an 
additive constant and A is unique. 

Now consider an arbitrary spherically symmetric (but 
not necessarily static) perturbation of the background 
metric. Let <fi(r, t) and X(r, t) denote, respectively, the 
perturbations of <E> and A. In other words, 2e _2 *</>(r, t) 
is the perturbation of the g metric component, and 
— 2e~ 2A A(r, t) is the perturbation of the g rr component. 
The gauge freedom in the perturbed metric is 6g ab — > 
Sg ab + C v (g^) ab , where (g ( - oy ) ab is the background met- 
ric and v a is an arbitrary vector field that generates dif- 
feomorphisms that preserve the metric form (|10|) . Since 
the only such nontrivial diffeomorphisms are t — > g(t), 
the only such nontrivial v a is v a — h(t)t a , where h(t) 
is a positive function and t a is the static Killing field of 
the background solution. Thus, we find that A is gauge 
invariant and <f> has the gauge freedom 

0(r,t)->#r,t) + /(t). (11) 

where / = dh/dt is an arbitrary function of t. Thus, 
our choice of the metric form (|10[) eliminates the gauge 
freedom in the metric perturbation variables except for 
the small residual freedom (|TT|> . 

An important consequence of the existence of this 
residual gauge freedom is that if A(r, t) and 4>{r,t) along 
with the appropriate perturbed matter variables solve the 
linearized equations of motion, then if we replace <p(r, t) 
by <f>(r,t) + f(t), we must still obtain a solution to the 
linearized equations of motion. This implies that <f> can 
appear in the equations of motion only in the form[iH| 
dcf)/dr. 



III. GAUGE FIXING 

As noted in the Introduction, generally covariant theo- 
ries are gauge theories: they possess "unphysical degrees 
of freedom" , corresponding to diffeomorphisms. In order 
to obtain deterministic equations of motion and obtain 
a variational principle of the sort we seek, we must elim- 
inate the gauge degrees of freedom. There is no known 
algorithm for doing this in general spacetimes. How- 
ever, in our work, we will be dealing with spherically 
symmetric perturbations of static, spherically symmet- 
ric spacetimes, and thus may restrict attention to space- 
times which are spherically symmetric but not necessarily 
static. For such spacetimes, it is always possible to put 
the metric in the following form: 

ds 2 = - exp(2$(r, t))At 2 + exp(2A(r, t))dr 2 + r 2 cK7 2 

^ (10) 
(See, e.g., [10( for the details.) Apart from the rotational 

isometries, the only diffeomorphisms that preserve the 

metric form (|10p are redefinitions of the time coordinate 

of the form t — > git), for an arbitrary monotonic function 

9- 



IV. SOLVING THE LINEARIZED 
CONSTRAINT EQUATIONS 

The general analysis of [9| shows that in any generally 
covariant theory, there will be a "constraint" on the phase 
space of the theory associated with the infinitesimal dif- 
feomorphism generated by an arbitrary vector field £ a . 
In other words, if £ a is used to define the notion of "time 
translations" , there will be a corresponding restriction on 
the phase space of the theory imposed by some subset of 
the equations of motion. If £ a is used to define the notion 
of "time translations" , then, as we shall see below, the 
resulting constraint equations will typically be of lower 
differential order in time than the other (so-called "evo- 
lution" ) equations of motion. Consequently, if we wish 
to get our equations of motion into the simple form |T]) , 
it normally will be necessary that we solve the constraint 
equations so that — after elimination of variables — the re- 
maining variables are unconstrained. Remarkably, we 
now shall show that in a wide class of diffeomorphism co- 
variant theories, for spherically symmetric perturbations 
of static spherically symmetric spacetimes, the linearized 
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constraint equations can be solved algebraically for the 
metric perturbation A in terms of the matter perturba- 
tion variables. Thus, for this wide class of theories, the 
constraints can be easily eliminated. 

Our analysis makes use of the fact proved in the Ap- 
pendix of 11] that for an arbitrary diffcomorphism co- 
variant theory, the Noether charge two-form can al- 
ways be defined so that that the Noether current (as 
defined in J7J), associated with an arbitrary vector field 
£ a takes the form 



CC a + dQ £ 



(121 



where, according to the analysis of [§], the equations 
£ a C a = are the constraints of the theory associated 
with the infinitesimal local symmetry £ a . (We will, in 
effect, rederive (p~2]) by our calculation leading to |28|) 
below.) Thus, if we perturb about a background solu- 
tion, we obviously obtain 



<5J ? = £ a 5C a + d(5Q € ) 



(13) 



On the other hand, we previously noted that 53^ satisfies 
@ . Now suppose that we are perturbing about a solution 
to the field equations, £ — 0. Suppose further that this 
background solution possesses a Killing field t a that is 
also a symmetry of all of the background matter fields, 
so that Ct^ — 0. Then, choosing £ a = t a , we find from 
CdD and © that 



t a SC a = d[t ■ 6 - SQ t ] 

Thus, the linearized constraint equations, t a SC a 
sociatcd with t a are equivalent to the equation 



where 



d/3 = 



P = t-G-6Q t . 



(14) 
0, as- 



(15) 



(16) 



The replacement of the equation t a 6C a = by the 
equation d/3 = need not, in general, result in any sim- 
plification of the equations, i.e., there may be as many 
or more independent components of (3 as t a SC a , and the 
equations d/3 = may be of as high or higher differential 
order as the equations t a SC a — 0. However, in the case 
of spherically symmetric perturbations of static spheri- 
cally symmetric solutions, the replacement of t a 5C a = 
by d/3 = always results in a major simplification. The 
reason for this simplification can be seen as follows. By 
spherical symmetry, in the coordinates introduced in the 
previous section, the three- form t a SC a must take the form 

t a 5C a = Hx(t,r)dt AdQ + H 2 (t,r)dr AdQ (17) 

where d£l = sin0 d9 Adip. Thus, the constraint equations 
give rise to two independent equations, namely Hi = 
and H2 = 0. By contrast, by spherical symmetry, the 
two-form (3 must take the form 



Thus, (3 has only one nonvanishing component, and the 
equation d/3 = then reduces simply to F = const. 
For the situations we shall consider here and in the sub- 
sequent paper, an "origin", r — 0, will be present in 
the spacetime. The two-form (3 is locally constructed 
from the background and perturbed dynamical fields 
(see (HE])), and thus must be smooth everywhere, includ- 
ing at r — 0. However, the spherical volume element 
dfl = sintfd^Ad^ is not smooth at r = 0. Consequently, 
we must have F — at r = 0, and since F is constant, 
we must have F = everywhere. In summary, we have 
shown that the two independent components of the lin- 
earized constraints take the form 



dF 



OF 

m,r) = w (19) 



and, thus, using the boundary conditions at r = 0, the 
two constraint equations Hi — and H2 = reduce to 
the single equation, F = 0. Furthermore, it clear from 
(fTi?]) that the equation F = will be of lower differen- 
tial order (by one) in the dynamical variables than the 
equations Hi = and Hi = 0. Thus, we obtain a ma- 
jor simplification by replacing the equations Hi = and 
H-2 = by the equation F = 0. 

It remains now only to get an explicit expression for 
F in an arbitrary diffcomorphism covariant theory. The 
calculation of is straightforward, and the expression 
for the Noether charge Q that satisfies (|T2"j) can be de- 
termined by following the procedures outlined in the Ap- 
pendix of [ll| ■ One may then obtain /3 — and, hence, F — 
from (fl6|) . However, instead of following this procedure, 
we shall give a very simple derivation of a general formula 
for the constraints t a C a . By linearizing this formula, one 
can then immediately determine Hi and H2 ■ The desired 
quantity F can then be determined by inspection from 
(O. 

As previously noted in Section [Til f° r an arbitrary vec- 
tor field £ a , we have 



dJ f = -£C e * ■ 



(20) 



Now, the right side of this equation is an four-form that 
depends linearly on £ a and contains no higher than first 
derivatives of £ a , i.e., it is of the form 



S£^ = {B a e + C\V a ^ b )e, 



(21) 



where e is the volume element associated with the met- 
ric and B a and C a b are locally constructed out of the 
dynamical fields. We can rewrite the right side as 



(B a C + C a b V a Z b )e = U a Ce + dV (22) 



where 



and 



U a = B a - W b C b 



(23) 



(3 = F(t,r) dn. 



(18) 



Vr 



<■<!( 



C a tacde 



(24) 
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Note that V does not depend upon derivatives of £ a . 
Thus, we have proven that for all £ a we have 



d(J 4 



CU a e. 



(25) 



However, the only way that this equation can hold for 
arbitrary £ a is if both sides are zero. Namely, if U a ^ at 
some point p, then we could find a smooth £ a of compact 
support such that U a £, a > everywhere and U a £ a > 
at p. However, the integral of the right side would then 
be positive, whereas the integral of the left side vanishes, 
thereby yielding a contradiction. Thus, we obtain 



and 



U a = B a - V b C b a = 



d( J £ - V) = . 



(26) 



(27) 



Equation (f26|) is a generalized version of the Bianchi iden- 
tity, applicable to an arbitrary diffeomorphism covariant 



theory, pos sibly containing matter fields. Equation (|27[ 
implies [12| that there exists an two-form Q' locally con- 
structed out of £ a and the dynamical fields such that 



J 5 = V + dQ' . 



(28) 



Comparing with (|12| and using the fact that V does not 
depend upon derivatives of £ a (the same argument as 
used below ([25)) ). we see that the constraints are given 
by 



ec a = V . 



(29) 



It is worth noting that, taking into account (|24|) and 
29| , the constraints take the form 



eC a =C-e, 



(30) 



where C a = C a b £, b , whereas the generalized Bianchi iden- 
tity (|26|) takes the form 



V&C a — B a 



(31) 



We will obtain explicit formulas for B a and C a b below. 
It should be emphasized that (f3"Tj) holds independently 
of whether the equations of motion are satisfied. 

We wish now to explicitly calculate V. For the pur- 
pose of this calculation, we assume that the Lagrangian 
C depends only on the inverse metric g ab and a sin- 
gle matter field A aia2 ~ an blb2 ,,, bm ; the generalization to 
the case of more than one matter field is straightfor- 
ward. We denote the gravitational equations of mo- 
tion (obtained by varying C with respect to g ab ) by 
(£G)ab£, and we denote the matter equations of motion 
(obtained by varying C with respect to A aia2 --- a " blb2 _ bm ) 
by (£M)a ia2 ...aJ lb2 - bm - Equation tfHD then takes the 
more explicit form 



dJ c = - {{£ G ) ab C ig ab 



+ (^M)a lC 



C^A a 



b 1 b 2 ...b„ 



(32) 



Using C^g ab = -(V a £ fc + V 6 £ Q ) and the standard formula 
for the Lie derivative of tensor field 

C 6 A a ^- a " blb2 ... bm = £ c V c A a ^- a " blb2 ... bm 

-^^••• c -°« blb2 ... fem v c r 

i 

+ J2A a ^- a « bl ... c ... bm V bi e, (33) 

i 

we can simply read off the formulas 

B c = -(£ M )a 1 ...a n bl - bm ^cA ai - a " bl ... bm (34) 

and 



C c d = 2g ca (£ G ) ad 

— ^ A ai '" an bl ,,,d...b m (^A/)ai...a„ &1 ' 



.c...b n 



J2 Aai - C - an b 1 ...b m (£M)a 1 ...d...aJ 1 - bm - (35) 



In the above equations, the summations run over all pos- 
sible substitutions of the indices c and d into the i th slot 
of A ai ° 2 - a » M2 ...6 m and (£M) ai a 2 ...a n blb2 - bm - We then 
see that 



U c = -2V a (£ G )ac ~ (£M) ai ...a n bl - bm VcA a ^ a " 



b, ...b„ 



-Y, V -A( £ M)a 1 ...c...a n bl --- bm A a ^-- a " bl ... bm ) 

i 

^V fc , (^ ai - a " 6l ... c ... 6m (£ A /) ai ...a„ bl - bm ) (36) 



and 



Vdef = Cede/ ^2^9 C (£ G ) ab 

+ J2^ Aai - an " 1 ...b m (£M)a 1 ...a n bl - 



c...b„ 



^r^ ai - C - Q "fc 1 ...b m (^/)a 1 ...a„ fcl - 6 "^ 



(37) 



Equation (|3"T|) is our desired formula for the constraints, 
which can be readily evaluated in cases of interest. 



V. ELIMINATING THE METRIC 
PERTURBATION VARIABLES 

We now wish to solve the linearized field equations 

(S£ G ) ab = 0, 6£ M = (38) 

for spherically symmetric perturbations of a static, spher- 
ically symmetric background solution. We shall show 
that, in a wide class of cases, it is possible to solve these 
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equations algebraically for the metric perturbation vari- 
ables, thereby reducing the problem to solving the lin- 
earized matter equations S£m — for the matter vari- 
ables alone. 

To begin, we note that from the form of C a b (see ([35])), 
together with the fact that the background equations of 
motion (£c)ab — 0, £m — are satisfied, it follows im- 
mediately that (|3"5|) is equivalent to 

5C ab = 0, S£ M = 0. (39) 

Next, we note that by spherical symmetry, there are only 
four independent components of 5C a b, which can be cho- 
sen to be the tt, tr, rr, and 99 components. However, we 
showed in the previous section that the tt and tr compo- 
nents of these equations can be replaced by a single equa- 
tion of lower differential order, namely F — 0. Further- 
more, the linearization of the generalized Bianchi identity 
([31) off of a background solution (where C a b = 0) yields 

V a {5C a b ) =8B b . (40) 

If we impose the linearized matter equations of motion, 
5£m = 0, then SBb = 0. This implies that if the tt, 
tr, and rr components of SC a b = are satisfied, and 
if the linearized matter equations of motion, S£m = 0, 
hold, then the 99 component of 8C a b = is automatically 
satisfied. Thus, the full set of linearized equations ([3"5|) 
may be replaced by 

F = , {5C) rr = , S£ M = (41) 

We now restrict consideration to theories in which the 
field equations are at most second order in derivatives of 
the metric, and such that the second derivatives of the 
metric appear only in the form of curvature. This is the 
case for general relativity with typical forms of matter, as 
well as for Bekenstein's TeVeS theory. It is also the case 
for f(R) gravity after a suitable redefinition of variables, 
since that theory can be recast as a scalar-tensor theory. 

We now analyze the possible dependence of SCu on 
A in order to determine the possible dependence of F 
on this quantity. By the Bianchi identity, the tt- and 
ir-components of the equations 5C ab = cannot con- 
tain second time derivatives of the metric perturbation 
quantities <fi and A, since otherwise the left side of (j40|) 
would contain third time derivatives of the metric pertur- 
bation quantities, which could not be canceled by terms 
on the right side, which, by assumption, contain at most 
second derivatives (see (|34|) ). (Recall that the gener- 
alized Bianchi identity (l3T|) holds independently of the 
field equations, so its linearization (|4"0")) off of a solution 
must hold if <j) and A are taken to be arbitrary func- 
tions of r and t.) Thus, in particular, the expression SCu 
cannot contain a term in d 2 X/dt 2 . On the other hand, 
for a diagonal metric like (fTDj) . for each coordinate x^, 
the Riemann tensor cannot contain a term of the form 
d 2 I 'dx^ 2 , since the formula for the Riemann tensor in- 
volves antisymmetrizations over components. Choosing 
x^ = r, we see that the Riemann tensor cannot depend 



on d 2 g rr I dr 2 . This implies that the expression SCu can- 
not contain a term of the form d 2 \/dr 2 . Finally, since 
the background spacetime is static and thus invariant un- 
der time reflection t — > —t, it follows that the expression 
SCu must not change sign under time reversal and thus 
cannot contain a term of the form dX/dt or d 2 X/drdt. 
Thus, we conclude that for theories in which the field 
equations are at most second order in derivatives of the 
metric, and are such that the second derivatives of the 
metric appear only in the form of curvature, the quantity 
A can appear in the expression SCu only in the form of 
A and dX/dr. However, as noted previously, the expres- 
sion F must be of lower differential order (by one) than 
SCu in all variables. Consequently, F can depend only 
algebraically on A. 

Next, we analyze the possible dependence of SCt r on 
<f> in order to determine the possible dependence of F 
on this quantity. As noted at the end of Section IIII| </> 
can appear in the linearized field equations only in the 
form of d(j)/dr and its derivatives. However, since SCtr 
must be odd under time reflection, (j> can appear in SCtr 
only in the form of d<fi/dt and its r-derivatives. Since, by 
assumption, the field equations contain no higher than 
second derivatives of the metric perturbation quantities, 
it follows that 4> can appear in 5C tr only in the form 
of d 2 4>/drdt. Suppose now that 8C tr contained a non- 
vanishing term proportional to d 2 4>/ drdt. Then, the r- 
component of the Bianchi identity (13 1 P would contain a 
term proportional to d 3 4>/ drdt 2 . The only other term in 
the Bianchi identity (|3ip that can contain third deriva- 
tives of the metric perturbation quantities is d(SC r r )/dr. 
In order to cancel the term proportional to d 3 (j>/drdt 2 , 
it is necessary that SC rr contain a term proportional to 
d 2 4>/dt 2 . However, this is impossible, since the field equa- 
tions can contain <f> only in the form of d(j)/dr and its 
derivatives. Consequently, SCtr cannot have any depen- 
dence whatsoever upon (f>. It then follows immediately 
that F cannot depend upon (f>. 

Putting together the results of the previous two para- 
graphs, we conclude that for theories in which the field 
equations are at most second order in derivatives of the 
metric, and are such that the second derivatives of the 
metric appear only in the form of curvature, we can solve 
the equation F — algebraically for X in terms of the 
matter variables. 

Next, we consider the equation (SC) rr — 0. For ex- 
actly the same reason as (5C) t t cannot contain second 
time derivatives of the metric perturbation quantities, 
it follows that (SC) rr cannot contain second derivatives 
of these quantities with respect to r, and thus cannot 
contain a term of the form d 2 (f>/dr 2 . By time reflec- 
tion symmetry, it also cannot contain a term of the form 
dcf>/dt or d 2 cf>/drdt. However, as already noted above, <f> 
can appear in the linearized field equations only in the 
form of dcf>/dr and its derivatives. Thus, wc also can- 
not have a term in [5C) rr of the form d 2 <t>/dt 2 (as we 
already noted above) or <fi. Consequently, (8C) rr can de- 
pend only algebraically on d<f>/dr. Thus, since we have 
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already eliminated A in terms of the matter variables for 
the class of theories considered here, we can solve the 
equation {5C) rr = algebraically for dtfi/dr in terms of 
the matter variables. 

In summary, we have seen that we can eliminate the 
metric perturbation variables algebraically by solving the 
equation F = for A and solving the equation (SC) rr = 
for d(f>/dr. We may then substitute these solutions into 
the linearized matter equations 



5£ 



M 







(42) 



to obtain a system of equations involving only the un- 
known matter variables. The perturbation problem has 
thus been reduced to solving these equations. 



VI. OBTAINING A VARIATIONAL PRINCIPLE 

Having fixed the gauge and completely eliminated the 
remaining metric perturbation variables by the proce- 
dures of the previous three sections, we may find that 
the system of equations (J42J) for the matter variables can 
be put in the general form ([1]). Indeed, we shall see in 
the next section and in a subsequent paper that this is 
the case in a wide class of theories. Nevertheless, even 
if the reduced equations take the form (fT]), our work is 
not completed because it is not straightforward to deter- 
mine if there exists an inner product that satisfies ([2]), 
and — even if such an inner product does exist — it is not 
straightforward to find it explicitly. This is particularly 
true if |T]) is a complicated system of equations. In the 
absence of an inner product that makes T a p self-adjoint, 
we will not be able to formulate a variational principle, 
and it will not be straightforward to analyze stability. 

We now shall show that the fact that the original the- 
ory was derived from a Lagrangian will enable us — at 
least in a very wide class of cases — to determine whether 
the desired inner product exists and to find it explicitly 
if it does. When the desired inner product does exist, 
we will thereby be able to immediately write down the 
variational principle ((3|) for determining stability. 

At first sight, it might appear that the fact that the 
original theory was derived from a Lagrangian would be 
of little use. It is true that one can obtain a Lagrangian 
for the linearized theory by expanding the Lagrangian 
of the exact theory to quadratic order about the back- 
ground solution. However, if one substitutes the gauge 
choice made in Section IIIII into this Lagrangian, the re- 
sulting object no longer functions as a Lagrangian. Sim- 
ilarly, even if one found a Lagrangian for the gauge-fixed 
theory, when one substitutes the solution for the metric 
perturbation variables found in Section |V] into this La- 
grangian, the resulting object would again fail to function 
as a Lagrangian. Thus, it is far from obvious how to ob- 
tain a Lagrangian for the reduced theory. 

Nevertheless, the fact that the original theory was de- 
rived from a Lagrangian leaves an important imprint on 
the reduced theory. As discussed in Section |TT1 in the 



original, exact theory, one can define a symplectic current 
three-form o;(^'; Si^f, 62^2), which is constructed out of a 
background solution ^ and two linearized perturbations, 
Si^k and 62^, such that u> depends linearly on 8i^> and 
62 Furthermore, the symplectic current satisfies the 
property that it is conserved, i.e., 



du; = 



(43) 



whenever Si*!/ and 62 satisfy the linearized equations 
of motion. Substitution of gauge choices for di^Sf and 
82^ and/or elimination of variables via some of the lin- 
earized field equations will not affect the conservation of 
u>. Thus, we automatically obtain a conserved symplectic 
current for the reduced theory. 

The conditions arising from the existence of a con- 
served symplectic current u; are most conveniently for- 
mulated in terms of the pullback, u>, of u; to the static 
hypersurfaces of the background solution, i.e., the hyper- 
surfaces orthogonal to the static Killing field t a . Since do; 
is a four-form in a four-dimensional space, the condition 
that do; = is equivalent to t ■ du; = 0. By a standard 
identity, we have 



t ■ du; = £ t u; — d(t • u;) 



(44) 



Thus, when the equations of motion hold, we have 

Ctu = d(t ■ uj) . (45) 

Pulling this equation back to the static hypersurfaces, we 
obtain 



£ t u; = d(t ■ ijj) . 



(46) 



In the theories we shall consider here and in a subse- 
quent paper, the pullback of the symplectic current for 
the reduced theory takes the form 



(47) 



where ip a denotes the dynamical variables for the reduced 
theory. Here, the three-form W Q/ 3 is constructed from 
the quantities appearing in the background solution, and 
thus is independent of t, i.e., £ t *W a p = 0. Thus, we 
obtain 



C t Q = W. 



a/3 



~df~df 







~di dt 



dt 2 



^ at 2 



(48) 



However, the linearized equations of motion (fTJ hold if 
and only if d 2 ip a /dt 2 = T a ^ . Thus, we find that the 
quantity 



W / 



'dipt <9i/>2 dtp% 0Vf 



dt dt 



dt dt 



(49) 
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must be an exact form on any static hypersurface, E t . 
Now, -0", dipf/dt, V>2) an d 9tp2 /9t are all freely specifi- 
able initial data at t = for the system of equations |TJ. 
Thus, if we choose , dipf /dt, and dip^/dt to be 
arbitrary smooth functions of compact support (or of suf- 
ficiently rapid decay), the integral of (|49|) over So must 
vanish. Inspecting the first two terms of (|49|) , we see that 
W Q ^ must be symmetric in its "field space indices", i.e., 

W Q/3 = W,g Q (50) 

and inspecting the last two terms of we see that 

for all ipi and 4>2 of compact support (or of sufficiently 
rapid decay), we have 

/ W Q ^T^ 2 7 " V> 2 Q ^V7) = 0. (51) 

Now, suppose that W Q/ 3 is positive definite in the sense 
that 

/ W a ^ a ^>0 (52) 

•/So 

for all ip a ^ 0. Then 

(V>i,V 2 )= / Wartfitf (53) 

defines an inner product. Equation (|5ip is then pre- 
cisely the statement that the operator T a p is symmetric 
in this inner product. We thereby obtain a variational 
principle [13], as desired. 

On the other hand, if W Q/ g fails to be positive defi- 
nite (or negative definite), then there docs not appear 
to be any reason to expect that an inner product exists 
that makes T a p self-adjoint [l^. In that case, we do not 
expect that there is a variational principle to determine 
stability, and one presumably must work directly with 
the equations of motion to analyze stability. 

Finally, it is worth noting that if ~W a p is positive def- 
inite, then given that T a p must satisfy ([ST) , it is easily 
verified that 

defines a Hamiltonian for the dynamics given by JTJ as- 
sociated with the symplectic form 

n(^,^)= / w(*;Vi,tfa)- (55) 

VII. ILLUSTRATION: CHANDRASEKHAR'S 
VARIATIONAL PRINCIPLE 

As a concrete example of the procedure outlined here, 
let us consider Einstein gravity minimally coupled to an 
isentropic perfect fluid. A variational principle for the 



radial oscillations of static, spherically symmetric stars 
was derived by Chandrasekhar [1, 0] by a direct analy- 
sis of the equations of motion. We now shall re-derive 
this variational principle in a much more systematic and 
direct way using the method described in this paper. 

In order to apply our method, it is essential that the 
equations of motion be derived from a Lagrangian. For 
general matter minimally coupled to g ab , the Lagrangian 
will be of the form 

C = -^Re + C mat [*,g ab } (56) 

where A denotes the collection of matter fields, with ten- 
sor indices suppressed. For the case of a perfect fluid, 
there have been there are several approaches that have 
been used to provide a Lagrangian formulation; see [l3j 
for an overview. We will use the "Lagrangian coordi- 
nate" method. In this formalism, we introduce an ab- 
stract three-dimensional manifold, M, of fluid worldlincs, 
equipped with a volume three- form N. The fluid is then 
described by a map \ '■ M — > M., which assigns to each 
x in the spacetime manifold M the fluid worldline that 
passes through x. By introducing coordinates, X A , on 
A4, where A = 1,2,3, we can represent x by the three 
scalar functions X A (x) on M, which we view as the in- 
dependent dynamical variables of the fluid. 

We define a three-form N a b c on spacetime — 
representing the "density of fluid worldlines" or the "par- 
ticle number density" — by 

N abc = N ABC {X)V a X A V b X B V c X c . (57) 

In terms of N abc we define the scalar particle number 
density v by 

v 2 - ^ Q6c jV afcc , (58) 
6 

and we define the fluid four-velocity U a by 

N abc = VtabcdU 11 , (59) 

where e a bcd is the spacetime volume four-form. The La- 
grangian four-form for the perfect fluid is then given by 

£mat = -Q{y)e, (60) 

where g(v) is an arbitrary function of the comoving par- 
ticle density v. The choice of the function g{y) corre- 
sponds to the choice of equation of state of the fluid (see 
(IB9|) below). 

We will need both the equations of motion and the 
formula for oj in our analysis, so we now proceed to 
derive these. As usual, the "gravitational part" of the 
Lagrangian, Cq = (l/167r)i?e, contributes (l/16n)G a b 
to the gravitational equations of motion and contributes 

m 

<4.v = S a bc d ef (S 2 g bc V d S ig e f - S ig bc V d S 2 g e f) , (61) 
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to the symplectic current, where 



d 

be e 



( ra cd ^ ad ^ ea 

/ = [OeO c9bf ~ 2# 9be9cf - b d c g e f 



1 



1 



2 S a eS d f9bc + -g ad g bc g ef ) (62) 

and we have defined u> a , the dual to u>, such that 

Wbcd — u a e abcd (63) 

with tubed defined by ([5]). 

To obtain the matter equations of motion and the con- 
tribution of the matter fields to to uj a , we vary £ ma t with 
respect to the dynamical fields X A . The variation of the 
number density three- form is given by (see [13]) 



SN abc = 3V [a (N ABC 6X A )V b X B V c] X c 



(64) 



(The antisymmetrization in this expression is, of course, 
over the spacetime indices only, not the fluid-space in- 
dices. Note that the fluid-space coordinates X A are sim- 
ply scalars as far as the spacetime derivative operator 
V a is concerned, and thus no metric variations, Sg ab , ap- 
pear in this expression. The antisymmetrization in this 
expression) Using ([58]) and 5v = 8(v 2 )/2v : we obtain 



6(-ge) = 



-£-V a (N ABC 6X A )V b X B V c X c N abc 
2v 



N ac dN b cd - gg ab I Sg 



2v 



,,ab 



(65) 



where g' denotes the derivative of g{v) with respect 
to v. The second term is brackets is the functional 
derivative of C mat with respect to g ab and thus is equal 
to —5 times the matter stress-energy tensor, T ab \ this 
quantity provides the contribution of £ ma t to the grav- 
itational field equations. Using (|59[) and the identity 
e abef e cdef = -^ c a S d \, we obtain 



N acd N b cd = 2v*{g ab + U a U b ), 



(66) 



and thus the stress-energy tensor can be rewritten in 
more familiar terms as 



T ab = g'vU a U b + {gv - g)g ab . 



(67) 



This can be recognized as the standard stress-energy for 
a perfect fluid 

T ab = + P)U a U b + Pg ab (68) 
under the identifications 

g —>■ p g'v — g — > P. (69) 

Integrating the first term in (f65|) by parts, we obtain 



V a (N ABC 5X A )\7 b X B V c X c N ab ' 



2;/ 



Va ( ~^-N ABC 5X A V b X B \7 c X c N abc 
2v 



+ N ABC 6X A X7 a ^X7 b X B X7 c X c N abc ) (70) 

From the first term on the right-hand side, we can read 
off the presymplectic current 9f nat : 



9mat - -fN ABC 5X A V b X B V c X c N abc 
2v 



(71) 



We can then take the antisymmetrized second variation 
°f $mat = • e, as in. |(5}, to obtain an expression for 
the symplectic current three-form w mat . By a straight- 
forward calculation, we obtain 

<at = [Sig bc K a bcA + 8 x X B L a AB 

+V b 5iX B M ab AB ] 8 2 X A -[1^2], (72) 

with w^ at defined analogously to ([6"3"} and the tensors 
K a bcA , L a AB , and M ab AB defined by 



K\ cA = ( ~ ^g bc N ADE V d X D V e X E N ade + 2N ABD V b X B V d X D N a c d + S a b N c de N ADE V d X D V e X E j 

+ ^3 (A ~ Q')N bde N c de N AFG V f X F V g X G N a ^, (73) 



L a AB = ~ (d B N ACD VcX c V d X D N acd + ZNacdV C X C V d X D d E N BFG V^X E V c X F V d ^X G ^ 

+ ^{e"v - g')N ACD V b X c V c X D N abc d E N BFG V e X E V f X F V g X G N e fv, (74) 
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and 

M ab AB = -|- (2N ABC W c X c N abc + 3N ACD W c X c W d X D N BEF g b ^W c X E W^X F ) 

+ -^{d'v - g')(N ACD V c X c V d X D N acd )(N BEF V e X E V f X F N b ^) . (75) 



(The antisymmetrizations in (|74p and (|75p are again over 
tensor indices only.) 

The second term in (1701) yields the matter equations 
of motion. Since the Lagrangian coordinates X A are 
scalars, we have S7[ a Vb]X A — 0, and since N abc is com- 
pletely antisymmetric, the vanishing of the second term 
for arbitrary SX A is equivalent to 

= V a (J-V b X B V c X°N abc ^ 

= ^\7 b X B X7 c X c X7 a {^N abc ^j . (76) 

This can be put in a more recognizable form by rewriting 
it in terms of the four-velocity U a 

- l -V b X B V c X c V [a (g'U d] ) e abcd . (77) 

This, in turn, is equivalent to U a V [„,([)' Ud\) — 0, which 
is just the relativistic Euler equation |13j . 

We now restrict our attention to spherically symmetric 
perturbations of static, spherically symmetric solutions. 
In order to obtain our desired variational principle, we 
must do the following: 

(i) Define our choice of variable(s) to describe the fluid 
perturbations; 

(ii) Obtain an expression for F and solve the equation 
F = algebraically for A; 

(iii) Write down the equation SC rr = and solve this 
equation algebraically for d<p/dr\ 

(iv) Substitute the solutions for A and d<j)/dr into 
the linearization of the matter equations of mo- 
tion (|76p . thereby rewriting this equation purely 
in terms of the perturbed fluid variable(s); 

(v) Determine if this equation takes the form dXJ) and, 
if so, read off the operator T"^; 

(vi) Evaluate the pullback of u> — or, equivalcntly, the 
time component of u a ; 

(vii) Determine if uj takes the form (|47p. and, if so, read 
off W a( 3 ; and 

(viii) Determine if W aj a defines a positive definite inner 
product and, if so, write down the variational prin- 
ciple ©. 



Although each of these steps may require some algebra, 
we emphasize that none of the steps require any inge- 
nuity, and we are guaranteed to succeed in obtaining a 
variational principle unless the matter equations of mo- 
tion obtained in step (iv) fail to take the form {TJ), the 
pullback of us fails to take the form (|47p , or W Q/ 3 fails to 
define a positive definite inner product. 

To describe the fluid perturbations, we choose the co- 
ordinates X A = {X R 1 X e , X®} on fluid space so that in 
the static background solution we have X R = r, X e = 9, 
and X^ = ip. Since we consider only spherically symmet- 
ric perturbations, we have SX — and 5X® = 0. Thus, 
the perturbation of the fluid is completely characterized 
by its radial "Lagrangian displacement" 

£(r,t) = 8X R (r,t) . (78) 

which describes the radial displacement of each fluid el- 
ement from its "equilibrium position" . 

Given our choice of X A , in order to be compatible 
with our assumption of spherical symmetry, the three- 
form N AB c on "fluid space" must take the form 

N Re *=q(X R )smX e , (79) 

for some function q. By (|58p . the number density, v, of 
the fluid is then given by 




The variation of this formula yields 

Sv-v(^-+(-— + — + ^\i-X\ (81) 
\dr \ v dr dr r J J 

This formula will enable us to express all terms in the lin- 
earized equations of motion that arise from the variation 
of v in terms of our chosen dynamical variable £. 

Our next step is to obtain an expression for F and 
solve the equation F = 0. Since, in this case, our matter 
fields X A are scalars, the terms in (|37p that depend on 
the matter equations of motion vanish, and we simply 
have 

V a bc — ^ d abct e {^G)de- (82) 

From this, we can conclude that 

V t6v = -2j—g(£ G ) tr g r \ V r g v = 2^=g~(£ G ) tt g tt ■ 

(83) 
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Since 5V = t a SC a , we can therefore identify the quantity 
Hi defined in (JT7J) as 



Hi(t,r) = -2r 2 e' s> - A (5£e) 



(84) 



Calculating the first-order equation of motion (6£a)trt 
we find that 



2<9A 



,2A 



From (|19p , we can immediately deduce that 

2 



F{t,r) =2r 2 e' p - A d< 



-2A 



A = 0. 



(85) 



(86) 



Thus, A can be eliminated in terms of the matter variable 
£by 



(87) 



The next step is to write down and solve the equation 
(SC) rr = for d(p/dr. By direct calculation, we obtain 

(SC) rr — (5£c)rr 

r \ or \ or r J I 



where 8v is given by (81). Thus, the solution to [SC) rr 
is[H 



dr \ dr dr J \ dr r 



~e 2 VV 



dr r v dr 



(89) 



Next, we need to write down the linearized matter 
equations of motion and substitute our solutions, (|87[) 
and (|89|) . for A and d<j)/dr into these equations to ob- 
tain equations written purely in terms of the perturbed 
matter variables. In our case, the linearization of (|77[) 
yields only one non-trivial equation for our single non- 
trivial perturbed matter variable £. By direct calcula- 
tion, we obtain 



„2A-2$ d 2 £ 
dt 2 



d$_ 

Or 



d_ 

()r 



d<f> 
dr 



q v 



dr 



1 dv 

v dr 



dA 
dr 



e-A 



= 0. 



(90) 



Substituting our solutions for A and d(j)/dr, we obtain 



allcl to t a , so t a N abc = vt a U d e abcd = 0. We obtain 



2A-2* 



d^_ 

dt 2 



-2*-A 



dr 



d<P 
dr 

3*+A 



d<S> 

dr 



dA 
dr 



II w ( 2 -* t-\ 



0. (91) 



Fortunately, this equation takes the form of (fTJ), and it 
is straightforward to read off the operator T from this 
equation. 

The next step is to evaluate the pullback of the sym- 
plectic current — or, equivalently, t a uj a — so as to obtain 
the inner product that will appear in the variational 
principle. As previously noted u> a consists of two pieces, 
u) a = cjg rav + Ci^t, given by (fBTjl and ([72]) respectively. 
By direct evaluation, we find that for spherically sym- 
metric perturbations of a static, spherically symmetric 
background, we have t a Wg rav = 0, i.e., there is no "grav- 
itational contribution" to the symplectic form[20]. In a 
static background, the expression for taU)^^ simplifies 
considerably, since the fluid four-velocity U a will be par- 



6 o 

2v 



N ABC V b X B V c X c S 2 X A [5ig ad N d 



be 



+3V [a (N DEF SiX D )V b X E V c] X F - [1 «-» 2]. (92) 

The first term in the square brackets also vanishes in 
the spherically symmetric case. Thus, we find that for 
spherically symmetric perturbations of static, spherically 
symmetric backgrounds, the symplectic form 



n 



d 3 xVX(< av +< lat )n a , (93) 



(where E is a static slice in the background spacetime, n a 
is its future-directed unit normal, and Vh is the volume 
element associated with the induced Riemannian metric 
h a b on E) is given by 

n = l J ^ x ^ n a^{N ABC 5 2 X A V b X B V c X c ) 

x (N DE F^ [a hX D ^ b X E V c] X F \ -[!<-> 2]. (94) 
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In writing the above, we have used the fact that there 
is only one non- vanishing component of 5X A , so any 
term proportional to SiX A S 2 X B (as opposed to terms 
depending on the derivatives of SX A ) will vanish under 
antisymmctrization. Writing this out in terms of our per- 
turbational variables yields the simple expression 

= tej drr 2 e 3A -V* (^6 - 6^) ■ 05) 
The inner product to be used in our variational principle 



can now be read off by comparing (|95|) with ([53]) . We 
obtain 

(6,6) = 4tt J dr r 2 e 3A -*g'vte 2 (96) 

Since q'v = p + P under the substitutions in ([69]) . this 
quadratic form will be positive for any fluid satisfying the 
null energy condition (i.e., p + P > 0.) 

Our variational principle is then of the form ([3]), with 
numerator 



(£,T£) = 127r J dr 



-r 2 e A " 



dr r J \ dr dr J 



'0 



(97) 



and denominator 

(£,£) = 12tt y drrV^V^ 2 (98) 

This is equivalent to the variational principle originally 
derived by Chandrasekhar [2l| 0, 0], but we have red- 
erived it here in a systematic way that is essentially "fool- 
proof" . 



IVII1 while the method can fail at various points (e.g., the 
inner product obtained from the symplectic form could 
fail to be positive), there is in principle never any "art" 
involved in applying this procedure to a given field the- 
ory. This method is therefore a potentially powerful tool 
for analyzing the viability alternative theories of gravity 
and other covariant field theories; in an upcoming work 
[l5|], we will apply this method to three alternative the- 
ories of gravity [J, H, Q . 



VIII. CONCLUSION 



We have presented a general procedure to analyse the 
stability of spherically symmetric perturbations about 
static spherically symmetric solutions of an arbitrary co- 
variant field theory. This procedure involves solving the 
linearized constraints and eliminating the metric pertur- 
bation variables algebraically. The symplectic form is 
then used to define an inner product, from which a vari- 
ational principle can be obtained. 

It is important to emphasize that our procedure is en- 
tirely prescriptive, as illustrated by the outline in Section 
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